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ABSTRACT 


Modern  optimal  control  methods  are  used  to  develop  a  three  controller 
system.  The  control  problem  Is  formulated  in  state  vector  form  using  full 
state  feedback.  The  control  is  implemented  on  a  modal  representation  of 
the  structural  equations  of  motion.  State  estimates  are  provided  by  a 
deterministic  observer.  A  computer  model  is  used  to  generate  controller 
and  estimator  gains  using  steady  state  optimal  regulator  theory.  Means  for 
extracting  structural  response  and  control  effort  time  histories  are 
developed  and  related  to  the  computer  model.  ^  ^ 

The  three  controller  system  is  used  to  control  the  CSDLI  model  sub¬ 
jected  to  physical  initial  conditions.  This  model  is  a  lumped  mass 
tetrahedral  truss  consisting  of  four  unit  masses  that  are  pin-connected 
by  massless  rods.  Two  versions  of  the  model  are  used:  The  nominal  ver¬ 
sion  as  designed  by  the  Draper  Lab  and  a  perturbed  version  designed  to 
minimize  the  deflection  due  to  a  static  load. 

Four  structure/controller  configurations  are  tested  via  computer 
simulations.  Compatibility  of  modal  groupings  and  the  use  or  non-use  of 
system  triangularization  are  the  factors  varied  to  form  the  controller 
configurations.  All  configurations  use  three  controllers  each  controlling 
four  of  the  twelve  modes.  Modal  input  to  the  program  is  supplied  by  a 
NASTRAN  analysis  of  each  structure.  Time  histories  of  structural  response 
and  expended  control  effort  are  computed  for  each  test  configuration  using 
a  series  of  structural  response  weightings  to  vary  the  control  and  estimator 
gains.  Performance  curves  are  formed  for  each  configuration  to  allow  com¬ 
parison  of  system  efficiencies.  The  effect  of  the  separate  combination  of 
modal  compatibility  and  use  or  non-use  of  system  triangularization  on  system 
efficiency  is  examined. '^iT^elative  ‘^controllability^  of  the  two  structures 
is  also  discussed. 


INVESTIGATION  OF  MULTIPLE  CONTROLLER 
SYSTEM  EFFICIENCY  AS  APPLIED  TO  THE  CSDLI  MODEL 

I.  INTRODUCTION 

Current  successes  and  projected  developments  in  the  space  transporta¬ 
tion  system  point  to  the  possibility  of  construction  in  space  of  structural 
systems  not  now  possible.  These  large  space  structures  (LSS)  may  be  too 
large  to  transport  and  too  flexible  to  withstand  earth's  higher  gravity 
environment.  The  projected  spans  of  these  systems  range  from  hundreds  to  ' 

thousands  of  feet.  Long,  lightweight  members  and  trusses  make  such  extreme¬ 
ly  large  structures  possible,  but  also  make  them  very  flexible.  This  added 

-  - 

flexibility  takes  the  form  of  an  increased  number  of  lower  frequency 
vibrational  modes.  Many  of  these  low  frequency  modes  will  fall  within  the 
bandwidth  of  the  control  system  and  thus  require  active  control. 

As  the  number  of  modes  that  must  be  included  in  the  control  model  -  i 

increases,  limits  in  on-line  computing  speed  and  off-line  computing  accu-  ] 

racy  force  the  use  of  multiple  controller  systems.  These  systems  are  j 

comprised  of  a  number  of  independent  controllers  each  controlling  a  sub-  -  -j 

set  of  the  controlled  modes.  The  number  of  modes  that  such  systems  can 
handle  seems  limited  only  by  the  expense  and  weight  of  the  extra  hardware 
and  by  the  ability  to  design  stable  systems.  Miller,  Calico  and  Miller, 
and  Aldridge  (Refs  1-3)  have  demonstrated  multiple  controller  design 
techniques  that  guarantee  overall  system  stability. 

As  larger,  more  flexible  structures  force  active  control  of  increasing 
numbers  of  modes,  the  increased  amount  of  energy  required  to  control  the 
system  becomes  a  factor  worthy  of  consideration.  Since  added  energy 
requirements  mean  additional  fuel,  electrical  sources,  or  solar  cell  area, 
the  problems  of  large  size  and  flexibility  could  be  further  compounded. 
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As  these  large  space  structures  become  reality,  so  will  the  need  to 
design  the  most  efficient  structure/control  system  configuration  possible. 
Current  structures  research  emphasizes  redesign  techniques  that  perturb  a 
baseline  design  to  yield  a  "more  controllable"  structure.  Since  much  of 
the  structural  geometry  and  placements  of  lumped  masses  on  a  structure  may 
be  predetermined,  most  of  these  methods  deal  with  techniques  for  deter¬ 
mining  stiffness  changes.  Design  methods  of  possible  use  include  those 
that  produce  structures  designed  for  separation  of  natural  frequ.  icies, 
minimum  deflection  under  a  static  load,  advantageous  modal  characteristics 
and  structural  approximation  of  the  optimal  regulator  control  solution. 
Because  structural  mass  is  always  a  consideration  in  evaluating  systems 
intended  for  use  in  space,  most  of  these  techniques  form  redesigns  holding 
the  total  mass  (or  equivalently,  the  total  stiffness)  constant.  It  is 
likely  that  the  most  advantageous  design  technique  for  a  particular  system 
will  depend  on  the  performance  criteria  used  to  evaluate  the  structural 
response. 

Proper  comparison  of  structural  designs,  or  in  general,  structural 
design  techniques,  requires  that  the  most  efficient  control  system  pos¬ 
sible  be  applied  to  each  structure  tested.  This  could  easily  require  a 
different  controller  configuration  for  each  structure  or  technique  com¬ 
pared.  Much  additional  study  is  needed  to  adequately  understand  the 
factors  involved  in  configuring  efficient  multiple  controller  systems. 
These  studies  will  require  a  computer  model  capable  of  forming  multiple 
controller  systems  and  simulating  the  response  of  the  control  system  as 
well  as  the  structural  response.  The  controller  configuration  must  be 
very  flexible,  and  the  simulations  should  be  capable  of  providing  any  data 
desired.  Aldridge  (Ref  3)  developed  a  computer  program  using  the  logical 


structure  presented  by  Miller  (Ref  1)  that  is  capable  of  forming  a  three 
or  four  controller  system.  The  program  can  also  decouple  the  controllers 
lO  form  a  decentralized  control  system  by  eliminating  the  combination  of 
observation  and  control  spillover  necessary  to  block  triangularize  the 
closed  loop  system  matrix. 

The  intent  of  this  thesis  is  to  investigate  the  effects  of  two  con- 
•.rj!  configuration  factors  on  multiple  controller  system  efficiency.  A 
•  4'^ied  version  of  Aldridge's  program  will  be  used  to  form  three- 
,ontrc;il*r  systems  for  the  controller  configurations  of  interest.  These 
.ontro!  i)$tems  will  be  applied  to  the  CSDLI  model  subject  to  a  set  of 
initial  conditions.  The  applicability  of  this  model  to  studies  of 
large  space  structure  control  problems  has  been  well  demonstrated.  Line 
of  sight  pointing  accuracy  is  used  as  the  criterion  of  system  performance. 
Control  system  performance  and  efficiency  will  be  evaluated  by  time 
histories  of  structural  response  and  total  expended  control  effort.  The 
controllers  designed  use  position  sensors  and  a  deterministic  observer  to 
estimate  modal  amplitudes  and  point  force  actuators  to  implement  the 
control . 

The  following  sections  detail  the  CSDLI  model,  the  computer  program 
wtth  the  modifications  used  to  form  the  controller  response,  and  the 
investigation  to  be  undertaken.  Following  the  discussion  of  the  investi¬ 
gation,  the  results,  conclusions  and  recommendations  are  presented. 


II.  TEST  MODEL 


The  demonstration  and  validation  of  structural  design  and  controller 
design  techniques  for  large  space  structures  is  a  difficult  task.  By 
nature,  large  space  structures  are  not  readily  categorized  by  a  few 
simple  examples.  This  is  largely  due  to  the  widely  varied  forms  of  large 
space  structures  and  the  large  number  of  structural  vibration  modes  which 
must  be  considered  in  modelling  them.  In  response  to  the  need  for  a 
simple  test  example,  the  Charles  Stack  Draper  Laboratory  (CSDL),  Inc.,  of 
Cambridge  Massachusetts,  developed  two  conceptual  structural  models.  The 
first  of  these  is  a  lumped  mass  tetrahedral  truss  referred  to  as  the 
CSDLI  model.  The  second  model,  CSDLII,  is  a  larger  order  generic  optical 
space  platform.  Only  the  tetrahedron  model  will  be  used  in  this  investi¬ 
gation.  The  following  presents  the  CSDLI  model  and  its  pertinent 
characteristics. 

CSDLI  Model 

The  CSDLI  model  was  selected  for  its  simplicity,  low  order,  and 
similarity  to  basic  large  space  structures  from  both  a  structures  and  a 
controls  point  of  view.  Because  this  model  has  only  four  nodes  which  are 
allowed  to  move  the  visualization  of  time  responses  in  physical  coordi¬ 
nates  is  quite  simple.  This  structure  provides  a  low  order  model  to 
which  control  systems  are  easily  applied.  Low  order  also  allows  all 
modes  to  be  both  modelled  and  actively  controlled.  This  simplifies 
verification  of  the  computer  model  and  reduces  uncertainties  in  the 
results  when  the  investigation  is  focused  on  structural  design  rather 
than  design  of  the  control  system.  Also,  the  response  characteristics 
of  this  model  are  very  similar  to  those  observed  in  large  space  structures. 


This  has  been  demonstrated  by  the  effective  application  of  decentralized 
control  systems  by  Miller  (Ref  1)  and  Aldridge  (Ref  3)  to  the  CSDLI  model. 

The  finite  element  model  of  the  structure  is  shown  in  Fig.  1.  The 
structure  consists  of  twelve  members  connected  by  pinned  joints  at  ten 
nodes.  This  allows  axial  forces  in  the  members,  but  no  bending  moments. 
The  members  are  considered  massless.  Equal  masses  of  one  unit  each  are 
lumped  at  the  four  vertices  of  the  tetrahedron  (nodes  one  through  four). 
Each  of  the  four  masses  has  three  translational  degrees  of  freedom  giving 
the  system  model  a  total  of  twelve  structural  modes.  Nodes  five  through 
ten  anchor  the  three  right  angled  bipods  which  support  the  tetrahedral 
truss  above  the  ground.  These  six  nodes  are  fixed  so  that  the  ground 
plane  forms  a  frame  of  reference  from  which  deflection  criteria  may  be 
established  and  measured.  The  coordinates  of  the  ten  nodes  are  listed 
in  Table  I. 

TABLE  I 

CSDLI  NODE  COORDINATES 


NODE 

X 

Y 

Z 

1 

0.0 

0.0 

10.165 

2 

-5.0 

-2.887 

2.0 

3 

5.0 

-2.887 

2.0 

4 

0.0 

5.774 

2.0 

5 

-6.0 

-1.155 

0.0 

6 

-4.0 

-4.619 

0.0 

7 

4.0 

-4.619 

0.0 

8 

6.0 

-1.155 

0.0 

9 

-2.0 

5.774 

0.0 

10 

2.0 

5.774 

0.0 

5 


The  origin  of  the  coordinate  frame  lies  in  the  ground  plane  beneath  node 
one.  The  Y-axis  lies  in  the  ground  plane  along  the  projection  on  the 
ground  of  the  member  joining  nodes  one  and  four  and  the  X-axis  forms  a 
right-handed  set. 

Six  pair  of  collocated  actuators  and  sensors  are  placed  on  the  six 
ground  connected  members  to  implement  the  control  system.  Appendix  A 
presents  the  actuator  direction  matrix,  D,  which,  in  the  case  of  col¬ 
located  actuators  and  sensors,  is  also  the  transpose  of  the  sensor 
direction  matrix,  C.  Also  included  in  Appendix  A,  are  the  frequencies 
and  eigenvectors  for  the  two  variations  of  the  CSDLI  model  that  are  used 
in  this  study.  This  data  was  obtained  by  NASTRAN  eigenstructure  analyses. 

This  study  uses  line  of  sight,  based  on  the  deflection  of  node  one 
from  the  Z-axis,  as  its  measure  of  controller  performance.  Any  motion  of 
node  one  in  the  Z  direction  is  considered  unimportant  to  performance.  The 
initial  conditions  listed  in  Table  II  are  applied  to  the  models  and  time 
histories  of  the  system  responses  are  computed. 

The  program  developed  by  Aldridge  (Ref  3)  was  modified  to  include 
output  of  the  integrated  total  expended  control  effort  in  the  time  history. 
The  integrated  control  effort  as  it  relates  to  the  reduction  in  line  of 
sight  deflection  can  then  be  used  to  evaluate  control  system  "efficiency". 
The  line  of  sight  deflection  and  system  efficiency  performance  parameters 
are  discussed  more  fully  in  later  chapters. 
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■  '  ■ '  •  ■ ■  ■  ’1 

• 

TABLE  II 

B 

CSDLI  INITIAL  CONDITIONS  IN  PHYSICAL 

COORDINATES 

-  - 

■ 

MODE 

DISPLACEMENT(q) 

VELOCITY (q) 

1 

-.00255 

-.00180 

1 

2 

-.00082 

-.00039 

— 

3 

.00394 

-.00175 

4 

-.0039 

-.00589 

■ 

1 

5 

-.00458 

-.02141 

-  J 

6 

.00352 

.00740 

1 

7 

.00031 

-.00327 

1 

R 

8 

-.00379 

-.00568 

-  - 

9 

-.00028 

-.00363 

10 

-.00058 

-.01546 

i 

11 

-.00499 

-.01586 

. 

12 

.00015 

-.00398 

III.  THE  COMPUTER  MODEL 


In  order  to  test  new  theories,  to  validate  structural  or  controller 
design  techniques,  or  to  test  actual  designs,  a  computer  simulation  is 
needed  that  can  serve  as  a  testbed.  The  AC0SS2  program  developed  by 
E.S.  Aldridge  (Ref  3)  forms  an  ideal  base  on  which  to  build  such  a 
testbed.  This  program  receives  input  of  structural  and  control  config¬ 
uration  data.  From  this,  it  forms  a  multiple  controller  system  capable 
of  handling  the  large  number  of  modes  required  to  control  a  large,  flexi¬ 
ble  structure.  The  program  then  decentralizes  or  decouples  the  multiple 
control  system  into  independent  systems.  The  control  system  designed  by 
this  program  is  output  in  the  form  of  controller  and  estimator  (or 
observer)  gains.  Eigenvalue  analyses  of  the  total  system  both  before  and 
after  the  decoupling  are  also  computed.  With  the  input  of  modal  initial 
conditions,  a  time  history  of  the  structural  response  may  be  computed. 

In  the  attempt  to  develop  a  testbed  program  from  this  base,  several 
qualities  appear  to  be  desirable:  1)  flexibility  in  the  inputs,  2)  mini¬ 
mum  outside  computation,  and  3)  inclusion  of  the  machinery  necessary  to 
generate  any  information  that  might  be  useful.  Improvement  in  these 
three  areas  is  the  motivation  behind  the  modifications  that  are  described 
in  this  chapter. 

This  chapter  consists  of  two  sections.  The  first  contains  a  summary 
of  the  development  of  multiple  controller  systems  as  presented  in  Ref  3. 
The  second  draws  on  the  development  of  the  first  in  order  to  describe  the 
modifications  that  have  been  made  in  the  AC0SS2  program. 
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The  Multiple  Controller  System 

As  presented  by  Calico  and  Miller  (Refs  1  and  2),  the  system  model 
for  any  structure  may  be  developed  from  the  vibrational  equations  of 
motion  in  the  form 

Mq  +  Eq  +  Kq  =  Du  { 1 ) 

where 

M  =  nxn  symmetric  mass  matrix 
E  =  nxn  symmetric  damping  matrix 
K  =  nxn  symmetric  stiffness  matrix 
q  =  nxl  physical  coordinate  vector 
u  =  mxl  control  input  vector 

By  introducing  the  nxn  modal  matrix  for  Eq  1,  such  that 

q  =  4>  *  n 

where  n  is  the  n-vector  of  modal  coordinates,  and  premultiplying  by  <l>^, 
Eq  1  becomes 

4>^M<frn  +  4>^E<l>n  +  <l>^K4>n  =  <t>^Du  (2) 

Due  to  the  properties  of  the  modal  matrix,  Eq  2  simplifies  to 


The  and  f  ;  being  the  natural  frequencies  and  damping  coefficients 
respectively  of  the  individual  modes.  Equation  3  can  now  be  written  in 
state  space  form  as 

X  =  Ax  +  Bu  (4) 


10 


where 


A 


I 


\  2 
-  CO 


I  \ 

1  -2r‘^ 

\  I  \ 


and 


(5) 


Because  the  complete  state  is  not  normally  available,  Eq  4  is  supple¬ 
mented  by  the  sensor  output  equation 


y 


C  q 


(6) 


where  the  p  and  v  subscripts  designate  position  and  velocity.  In  terms  of 
the  state  vector  X  this  is  expressed  as 


y  =  Cx 


(7) 


where 

c  =  [Cp  i  C,  ]  (8) 

when  only  position  sensors  are  used,  Eq  8  becomes 


Further  if  collocated  actuators  and  sensors  are  used 

y  =  [d''’4;o]  =  0]  (10) 

The  program  reads  the  C  matrix  in  transposed  form  to  take  advantage  of 
this  simplicity,  but  does  not  require  that  C  equal  D  transported.  Equa¬ 
tions  4  and  9  form  the  state  space  model  that  is  used  to  design  the 
control  system. 

This  study  will  use  the  three  controller  program  configuration  with 
no  modes  left  uncontrolled  (no  residual  modes).  The  following  develop¬ 
ment  of  a  multiple  controller  system,  therefore,  assumes  this  configura¬ 
tion.  For  a  four  controller  system  or  three  controllers  with  residuals, 
the  reader  is  referred  to  Ref  3. 

For  three  controllers,  each  controlling  a  subset  n-j  of  the  con¬ 
trolled  modes,  the  state  vector  x  is  simply  represented  as 

X  =  (  X^ ,  X2,  X3  ]  (11 ) 

The  x^  terms  represent  2n^ -vectors  that  are  the  state  vectors  of  the 
three  controllers.  In  the  form  of  Eqs  4  and  9,  the  state  equation 
models  of  the  individual  controllers  are 

X.  =  A.x.  +  B.u  (12) 

and 

3 

y  =  s  C.x.  (13) 

i=l  ^  ^ 
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where  the  and  have  the  form  of  Eq  5  and  the  have  the  form 

of  Eq  9.  Now,  the  desired  state  feedback  control  has  the  form 

3 

u  =  E  G.x.  (14) 

i=l  ^  ^ 

where  are  the  control  gain  matrices. 

As  mentioned  earlier,  the  entire  state  vector  is  not  normally 
available,  so  an  estimator  is  needed  to  estimate  x  from  the  observations 
y.  This  estimator  takes  the  form 

X.  =  A.  X.  +  B.ii  +  K^.{y-y^)  (15) 

where 


A 


(16) 


and  -A  represents  an  estimated  vector.  The  K^.  are  the  estimator  gains  and 
are  chosen  so  that  the  state  estimate  error 


(17) 


is  asymptoticly  stable.  Now,  the  control  vector  is  redefined  in  terms  of 
the  estimated  state  by 

3 

u  =  E  G.  X.  (18) 

i=l  ^  ^ 


where 


A 


X.  .Ei 


(19) 


The  controller  and  estimator  gains  are  computed  using  linear  optimal 
regulator  theory.  Regulator  theory  is  developed  very  completely  in 


Kwakernaak  (Ref  4)  and  is  applied  specifically  to  our  case  in  (Ref  3). 

This  development  needs  only  the  fact  that  gains  found  by  regulator  theory 

minimize  the  quadratic  cost  function  (or  performance  index)  J  defined  by 

00 

J  =  1/2  J  (xTq^x^  +  u^R^.u)dt  (20) 


where 


Q  is  an  nxn  positive  semidefinite  response  weighting  matrix 
R  is  an  mxm  positive  definite  control  weighting  matrix  subject  to 
Eq  12.  In  the  AC0SS2  program,  the  R  matrices  are  all  the  identity  and 
the  response  weighting  is  restricted  to  being  diagonal  with  the  diagonal 
elements  input  by  mode. 

Once  the  control  and  estimator  gains  are  computed,  the  program  forms 
an  augmented  state  equation  that  may  be  integrated  to  form  the  time 
history  of  the  state  vector.  This  is  done  by  forming  an  augmented  state 


vector  z  given  by 


;  -  fjT  pT  -T  pT  -T  -i 

Z  -  I ,  ^2 j  Xj,  J 


-T  ,  T 


where 


'i  = 


Expressing  Eq  12  in  terms  of  z  and  forming  an  i;  Equation 
Eq  12  in  terms  of  state  z  elements  are 


+  6,G,i,  4  E  B,Gj(ij4ej)  t  =  1,2,3 

J  • 
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Differentiating  Eq  17  gives 


®i 


Xi  -  X. 


(24) 


Then  combining  Eqs  15,  16,  19,  and  23  yields 


i,  =  (Ai-K,C.)i,  *  I  K,C,Sj 

J  ■ 


3 

E 

i=i 

j+i 


i=l,2,3 


(25) 


In  matrix  form,  the  closed-loop  state  equation  in  terms  of  z  is 


z  = 


il+BiGi 

BiGi 

®iS 

B1G2 

B1G3 

B1G3 

0 

A^-KiCi 

K1C2 

0 

K1C3 

0 

B^Gi 

B^G, 

A2+B2G2 

^2^2 

B2G3 

^2^ 

K2C1 

0 

0 

A2-K2C2 

K2C3 

0 

¥1 

B3G1 

^3^2 

^3*^2 

^3^3 

^3^1 

0 

K3C2 

0 

0 

^3~^3^ 

(26) 


The  solution  to  this  differential  equation  is 
z  =  EXP[At]  *  i(o) 


(27) 


where  A  is  the  matrix  of  Eq  26.  This  is  the  solution  that  is  propagated 
to  form  the  response  time  history.  Note  that  knowledge  of  the  state  z 
is  not  a  very  useful  output.  The  state  z  is  a  vector  of  modal  coordinates 
reordered  into  three  controllers,  each  controller  state  being  augmented 
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with  an  error  term.  The  method  used  by  the  program  to  extract  useful 
information  is  discussed  in  the  following  section  of  this  chapter. 

Decentralization  or  decoupling  of  the  controllers  is  accomplished 
by  lower  block  triangularization  of  the  matrix  in  Eq  26.  This  is  done 
by  insuring  that  the  proper  combination  of  control  and  observation 
spillover  is  eliminated  by  the  decoupling  transformations  developed  in 
Ref  3.  As  discussed  in  the  Investigation  and  Results  section,  this 
transformation  is  not  always  necessary  and  it  may  reduce  the  efficiency 
of  the  control  system. 

Program  Modifications 

As  mentioned  earlier  in  this  chapter,  three  goals  of  this  modifica¬ 
tion  were  flexibility  in  the  input,  minimum  outside  computation,  and  the 
ability  to  generate  any  needed  output.  The  changes  made  in  the  program 
in  the  course  of  this  study  reflect  improvements  in  each  of  these  areas. 
Further  study,  however,  may  well  reveal  other  advantageous  revisions. 
Certainly,  the  unique  needs  of  particular  tests  will  require  small 
changes. 

The  following  presents  the  modifications  made  thus  far.  They  are 
discussed  in  the  order  that  the  program  encounters  them.  Because  only 
the  modifications  are  discussed,  ties  to  the  previous  section  will  be 
few,  but  equation  numbers  are  given  where  applicable.  The  variable  and 
matrix  names  used  by  the  program  are  either  used  in  the  discussion  or 
listed  in  parentheses.  The  FORTRAN  code  of  these  modifications,  as  well 
as  a  program  flowchart  and  a  sample  data  file  are  included  as  Appendix  B. 
Where  possible,  a  few  lines  before  and/or  after  the  change  are  listed  to 
show  location  in  the  program. 
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The  first  change  alters  the  way  in  which  the  structural  and  control 
hardware  information  are  input.  Because  sensor  and  actuator  location 
are  logical  test  variables,  direct  inputs  of  4^D  (Eq  2)  (PHIA)  and  C4 
(Eq  9)  (PHIS)  were  changed  to  allow  separate  inputs  of  4  (PHI),  D  (DMAT), 
and  C  (CMAT).  Matrices  PHIA  and  PHIS  are  then  computed  internal  to  the 
program.  Now  the  only  outside  computation  required  is  a  NASTkAN  analysis 
of  the  structure  and  the  assemblage  of  the  C  and  D  orientation  matrices. 
Here  a  new  subroutine,  TANGL,  computes  the  angles  between  the  rows  of 
PHIA.  As  mentioned  by  Miller  (Ref  1),  the  choice  of  modal  groupings 
within  individual  controllers  is  an  important  consideration.  According 
to  Miller,  the  angles  between  the  control  forces  applied  to  the  modes 
within  a  single  controller  should  be  kept  as  small  as  possible.  This 
point  is  the  subject  of  a  quantitative  investigation  later  in  this 
study. 

Another  change  in  the  input  is  the  addition  of  an  option  to  input 
the  initial  conditions  in  either  modal  coordinates  (INIT)  or  physical 
coordinates  (GNIT).  This  option  is  presently  assigned  by  initialization 
of  the  variable  INCOND  within  the  source  program.  For  the  purpose  of 
this  study,  equal  initial  conditions  in  physical  variables  were  required 
for  the  fair  comparison  of  two  structural  designs.  If  modal  input  is 
chosen,  the  physical  coordinates  are  computed  and  output  for  comparison. 
The  program  uses  modal  coordinates  to  form  the  controllers  and  the  time 
history  without  regard  to  the  choice  of  input. 

The  line  of  sight  matrix  (CLOS)  that  extracts  the  performance 
parameters  for  the  time  history  is  input  as  in  AC0SS2.  However,  it 
should  be  noted  here  that  its  use  has  been  altered  significantly.  This 
will  be  discussed  in  detail  later.  Input  of  the  diagonal  values  of  the 


weighting  matrix  (Eq  20)  (BB)  changed,  however,  these  values  now  set 
only  the  relative  weighting  between  the  modes.  The  magnitudes  of  these 
values  are  then  scaled  by  the  value  of  SCALEB  which  is  read  from  data 
during  the  program  initializations. 

The  only  change  in  the  section  that  forms  the  individual  controller 
matrices  (Eqs  12  &  13)  (A1-A3,  B1-B3,  C1-C3)  is  the  addition  of  the 
formation  of  PHI1-PHI3.  These  matrices  consist  of  columns  containing 
the  eigenvectors  of  the  individual  controllers.  They  are  generated  by 
subroutine  FORMPHI  and  are  used  later  in  the  program  to  form  the 
reordered  modal  matrix  (PHINEW). 

The  remainder  of  the  changes  deal  with  the  machinery  used  to  extract 
useful  time  histories  from  the  augmented  state  vector,  Z  (Eq  27).  These 
changes  are  contained  in  the  area  of  the  main  program  just  before  sub¬ 
routines  YHAT  and  TIME  and  in  the  subroutines  themselves.  The  program 
stores  the  matrix  of  Eq  26  in  array  MAJM  so  that  the  solution  to  Eq  27  is 

i  =  exp  [MAJM*t]  *  z(o)  (28) 

where 

z  —  ^n*j  ,n-j  je*!  ,e*|n2n2,e2*62,n2n2,e0,e2^ 

and  z(0)  has  been  formed  from  the  modal  initial  condition  matrix  (INIT) 
by  FORMX0.  The  state,  z(0),  is  then  propagated  by  subroutine  time  to 
arrive  at  z(t).  However,  as  previously  mentioned,  z(t)  is  a  combination 
of  modal  coordinates  and  their  rates  as  well  as  observation  error  and 
error  rate  terms.  The  program  uses  a  series  of  matrices  to  extract  the 
physical  coordinates,  q(t),  from  the  full  state,  z(t).  To  obtain  a 
reordered  set  of  modal  coordinates,  YHAT  forms  the  matrix  PETA  such  that 
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PETA  *  z  = 


(29) 


where  is  the  modal  coordinate  vector  for  the  i^”  controller.  The  main 
program  has  formed  the  matrix  PHINEW  of  reordered  mode  shapes  so  that 


q  =  4n  =  PHINEW  *  n^ 


where  q  is  the  physical  coordinate  vector  (Eq  1)  in  its  original  order. 
This  is  true  because  the  columns  of  PHI  were  reordered  to  form  PHINEW  in 
the  same  manner  that  the  rows  of  n  where  reordered  to  form 

The  structural  time  response  of  interest  to  the  designer  may  be 
individual  physical  coordinates,  a  linear  coordination  of  physical  coordi¬ 
nates,  or  a  nonlinear  function  of  coordinates.  The  matrix  CLOS  is  used 
by  the  program  to  form  desired  response  vector,  XLOS,  from  the  physical 
state  so  that 


XLOS  =  CLOS  *  q 


If  a  nonlinear  function  is  the  desired  performance  parameter,  CLOS  may  be 
used  to  extract  the  states  that  are  needed,  and  the  function  may  then  be 
computed.  The  present  study  uses  this  method  to  form  the  deflection 
parameter  which  is  a  nonlinear  function  of  the  x  and  y  coordinates  of 
node  one.  Combining  Eqs  29,  30,  and  31  yields 


XCLOS  =  CLOS  *  PHINEW  *  PETA  *  z 
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This  use  of  CLOS  differs  from  the  AC0SS2  program  which  uses  a  matrix 
V  in  place  of  PHINEW.  The  matrix  V  transforms  [  n|,  n^.n^  ]  into  the 
vector  y  of  sensor  output.  It  was  found  that  this  method  makes  the  for¬ 
mation  of  CLOS  unnecessarily  complicated  if,  indeed,  a  physical  response 
parameter  is  the  goal.  The  machinery  for  forming  V  was  not  removed,  how¬ 
ever,  so  that  the  sensor  output  time  history  can  be  formed  with  reasonable 
effort. 

Subroutine  THAT  multiples  CLOS  *  PHINEW  and  calls  the  product  CV,  forms 
PETA,  and  then  performs  CV*PETA  calling  the  new  matrix  XHAT.  Equation  32 
may  now  be  written 

XLOS  =  XHAT^i  (33_ 

which  is  the  form  used  by  subroutine  TIME  to  form  the  time  history. 

A  procedure  similar  to  the  formation  of  XLOS  is  used  to  form  the 
control  vector,  u  (UBAR).  Since 

3 

G  =  I  G.x.  (34) 

i=l  ’  ’ 

The  estimated  state,  x^  must  be  formed.  From  Eq  19 

,  "i  *  ^ 

fi,  *  i, 

To  extract  the  x.,  subroutine  THAT  forms  PXHT  so  that 


To  form  and  sum  the  products  ,  the  main  program  forms  a  partitioned 
matrix  of  controller  gains 

GPART  =  [G^l  Gg  ;  G3]  (37) 

Now  u  is  formed  easily  by 

3 

=  z  G.x.  (38) 

i=l  ^  ^ 

Combining  Eqs  36  and  38  yields 

u  =  UHAT  *  z  (39) 

where 

UHAT  =  GPART*PXHT  (40) 

With  UHAT  input  into  subroutine  TIME,  u(t)  is  readily  formed  from  the 
augmented  state,  i(t).  The  present  study,  however,  uses  the  integrated 
or  total  expended  control  force  as  the  parameter  of  interest.  To  perform 
the  integration  of  the  control  force,  a  trapezoidal  integration  scheme  is 
included  in  subroutine  TIME.  Direct  output  of  u(t)  is  possible  without 
removal  of  the  integration  section  by  simply  altering  the  WRITE  statement. 

With  this  machinery  in  place,  it  is  hoped  that  the  minor  modifications 
needed  to  suit  the  needs  of  a  particular  test  may  be  accomplished  with 
relative  ease.  At  this  time,  however,  the  warning  contained  in  the  recom¬ 
mendation  section  should  be  considered  the  probable  causes  of  any 
unexplainable  major  problems. 


IV.  INVESTIGATION 


Control  of  large  space  structures  is  a  complex  task  that  must  be 
carefully  dissected  to  be  fully  understood.  Complete  knowledge  of  the 
problem  and  the  possible  paths  to  the  "best  solution"  requires  a  thorough 
knowledge  of  the  individual  factors  involved  and  of  the  sensitivity  of 
each  factor  to  changes  in  the  others.  The  studies  done  by  Miller  and 
Aldridge  show  that  multiple  controller  systems  are  both  needed  and  viable 
as  solutions  to  the  space  structure  control  problem.  This  investigation 
has  two  thrusts;  first,  to  quantitatively  test  two  qualitative  insights 
mentioned  in  the  previous  studies.  These  deal  with  the  grouping  of  modes 
into  controllers  by  the  angles  between  the  vectors  of  modal  amplitudes 
(between  the  rows  of  4^D)  and  the  general  tendency  of  the  decoupling 
transformation  to  reduce  system  stability.  These  areas  are  studied 
individually  and  in  relation  to  one  another.  The  second  thrust  is  to 
demonstrate  a  method  of  comparing  the  efficiency  of  one  structure/ 
controller  design  with  that  of  another  design. 

Test  Cases 

Four  structure/controller  configurations  are  tested.  Three  are 
comprised  of  control  configuration  changes  on  the  nominal  structure  and 
the  fourth  uses  the  deflection  design  to  demonstrate  the  considerations 
necessary  to  fairly  compare  two  structural  designs. 

Several  factors  of  the  control  configuration  remain  unchanged  from 
case  to  case.  The  program  is  used  only  in  the  three  controller  configura¬ 
tion  with  four  of  the  twelve  structural  modes  in  each  controller.  A 
structural  damping  ratio  of  1/2%  (.005)  is  applied  to  each  mode.  The  six 
pair  of  collocated  actuators  and  sensors  are  located  on  the  six  members 


which  support  the  tetrahedral  truss.  The  modal  initial  conditions  used 
by  Miller  are  applied  to  the  cases  involving  the  nominal  structure.  The 
program  output  for  the  nominal  structure  includes  initial  conditions  in 
physical  coordinates.  These  values  are  used  for  the  case  using  the  fre¬ 
quency  designed  structure  so  that  the  physical  initial  conditions  are 
constant  over  the  four  cases. 

The  response  weighting  matrix,  Q,  is  the  varied  parameter  within  the 
tests,  but  the  values  used  are  the  same  for  the  four  tests.  Q  is  input  as 
the  identity  and  scaled  by  the  value  of  SCALEB  in  the  program.  Since 
specific  response  requirements  were  not  the  goal,  the  range  of  SCALEB  used 
in  the  tests  was  chosen  solely  for  clarity  of  presentation.  By  a  rough 
search  of  a  few  orders  of  magnitude,  a  suitable  range  was  found  to  be  0.1 
through  1.0.  To  obtain  clear  separation  between  runs,  the  values  used  were 
lE-10,  0.04,  0.1,  0.2,  0.5,  0.8,  and  1.0.  The  first  two  values  were  added 
to  complete  information  needed  to  form  the  plots  in  the  next  chapter. 

The  three  nominal  design  cases  are  separated  by  the  manner  in  which  the 
modes  are  assigned  to  the  controllers  and  by  whether  or  not  the  decoupling 
transformation  was  included  in  the  controller  design.  A  modal  grouping  that 
minimizes  the  angles  between  the  modes  within  a  controller  is  referred  to 
as  favorable.  An  unfavorable  grouping  assigns  orthogonal  modes  to  the  same 
controller  and  modes  with  a  small  relative  angle  to  separate  controllers. 
Certainly,  there  are  varying  degrees  of  unfavorable  groupings.  This  study, 
however,  uses  only  the  worst  available  grouping  to  represent  unfavorable. 

The  transformation  that  block  triangul arizes  the  system  matrix  has  been 
shown  to  lessen  the  damping  originally  achieved  by  the  controllers.  While 
the  transformation  guarantees  overall  system  stability,  it  is  quite  possible 
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that  some  control  configurations  are  stable  before  triangular! zed.  The 
configurations  that  were  tested  are  now  presented  case  by  case. 

Case  I.  The  nominal  structure  modes  are  grouped  favorably,  and 
the  decoupling  transformation  is  used.  Table  III  presents  the  relative 
angles  for  the  nominal  structure  modes  as  computed  by  the  program.  The 
most  favorable  grouping  is  given  by 
Controller  1 :  2, 3, 8, 9 

Controller  2:  5,6,11 ,12 

Controller  3:  1 ,4,7,10 

Case  II.  The  nominal  structure  modes  retain  the  favorable  grouping 
of  Case  I,  but  the  decoupling  transformation  is  omitted. 

Case  III.  The  nominal  structure  modes  are  grouped  unfavorably.  From 
Table  III,  one  of  the  most  unfavorable  groupings  is 
Controller  1 :  2, 3, 4, 6 

Controller  2:  1 ,5,7,9 

Controller  3:  8,10,11,12 

Of  the  six  modal  angles  in  each  controller.  Controller  1  has  five  that 
are  essentially  90°,  Controller  2  has  four  that  are  orthogonal,  and  Con¬ 
troller  3  has  five  angles  greater  than  87°.  The  decoupling  transforma¬ 
tion  is  included  in  the  controller  design  for  this  case. 

Case  IV.  The  deflection  designed  structure  modes  are  grouped 
favorably,  and  the  decoupling  transformation  is  performed.  This  case  is 
intended  to  represent  the  best  possible  control  configuration  for  the 
deflection  design.  The  reasons  for  its  choice  are  discussed  in  the  next 
section.  Table  IV  presents  the  angles  between  the  modes  of  this  structure. 
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TABLE  III 

RELATIVE  ANGLES*  BETWEEN  NOMINAL  DESIGN  MODES 


1 

o 

o 

89.99 

-90.00 

89.99 

0.00 

64.38 

®  3 

-90.00 

64.38 

0.00 

4 

66.58 

89.99 

-90.00 

5 

90.00 

-83.28 

85.33 

89.99 

89.99 

90.00 

■  7 

33.16 

90.00 

90.00 

8 

90.00 

-50.03 

68.61 

9 

89.99 

32.74 

34.86 

10 

-51.61 

89.99 

89.97 

^  11 

-90.00 

-85.04 

-61.81 

12 

89.99 

71.78 

-89.56 

p.  MODE  1 

1  7 

8 

9 

66.58 

89.99 

-90.00 

0.00 

-90.00 

90.00 

-80.25 

89.99 

89.98 
-14.98 
-90.00 

89.99 


90.00 
-83.28 
85.33 
-90.00 
0.00 
90.00 
90.00 
-84.41 
-77.88 
-89.99 
-23.55 
-n  .54 


9.99 

9.99 


90.00 

90.00 

90.00 

0.00 

-90.00 

-90.00 

90.00 

-90.00 

-90.00 

-90.00 


The  most  favorable  grouping,  though  not  so  obvious  as  for  the  nominal 
design,  appears  to  be 

Controller  1 :  3,4,5,11 

Controller  2:  1 ,2,6,12 

Controller  3:  7,8,9,10 

The  modified  program  is  run  seven  times  on  each  test  case  varying  only 
the  scale  of  the  response  weighting  matrix.  The  output  of  interest  is 
the  time  history  generated  by  application  of  the  given  initial  conditions 
for  each  case  and  weighting  matrix.  The  parameters  included  in  each  time 
history  are  the  x  and  y  coordinates  of  Node  1,  the  radius  of  deflection 
for  node  1,  the  integrated  effort  of  each  actuator,  and  the  total 
integrated  effort  for  all  actuators.  The  following  section  presents  and 
discusses  the  results  obtained  from  these  computer  runs. 
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VI.  RESULTS 


Six  time  response  simulations  were  performed  on  each  of  the  four  test 
cases  described  in  the  previous  section.  Different  controller  and  estimator 
gains  were  designed  for  each  simulation  within  a  test  case  by  varying  the 
response  weighting  matrix,  Q.  In  this  chapter,  the  factor  used  to  scale  Q 
will  be  referred  to  as  so  that  ctq  is  the  varied  parameter  within  a 
test  case.  In  addition,  uncontrolled  responses  to  the  initial  conditions 
were  approximated  for  both  structural  designs  by  setting  ag  equal  to 
lE-10.  This  value  of  aQ  resulted  in  zero  gains  for  all  controllers  and 
zero  expended  control  effort.  The  data  obtained  from  these  simulations 
is  presented  in  the  form  of  three  plots  for  each  test  case.  While  the 
plots  contain  axes  labelled  in  terms  of  force  and  distance,  no  units  are 
given.  This  is  because  the  units  are  of  no  consequence  to  this  study  as 
long  as  consistent  mass,  distance,  force,  and  stiffness  units  are  used 
so  that  the  natural  frequencies  are  in  hertz. 

The  first  plot  contains  the  response  envelopes  for  the  values  of 
that  were  used.  As  mentioned  earlier,  the  x-  and  y-  line  of  sight 
parameters  used  in  previous  studies  on  this  model  have  been  used  to  form 
a  single  response  parameter,  the  deflection  node  one  from  the  z-axis. 

The  use  of  a  single  response  parameter  simplifies  comparison  of  responses. 
Also,  the  radius  of  deflection  time  history  can  be  described  by  a  single 
response  envelope  because  it  is  always  a  positive  quantity.  Figure  2 
demonstrates  the  formation  of  a  typical  response  envelope.  It  should  be 
noted  that  the  determination  of  the  response  envelope  from  the  complete 
time  history  is  somewhat  subjective.  The  time  step  used  (0.01  sec)  was 
sufficiently  small  to  show  the  oscillations  completely.  Still,  the  peaks 


of  individual  oscillations  were  not  always  well  defined.  The  height  of 
peaks  located  in  the  center  of  a  time  interval  was  impossible  to  determine 
accurately.  This  problem  is  illustrated  clearly  on  Figure  2  at  about  0.3 
seconds.  When  the  height  of  an  oscillation  was  very  uncertain,  the  peak 
was  not  plotted.  Each  point  plotted  is  subject  to  some  levels  of  error 
due  to  the  discrete  nature  of  the  time  history.  Because  these  errors 
always  cause  a  lowering  of  the  peaks  from  the  actual,  the  response  enve¬ 
lope  was  not  allowed  to  fall  below  any  point  of  the  time  history.  Scatter 
in  the  data  is  due  to  these  errors. 

This  problem  was  magnified  on  the  frequency  design  to  the  extent  that 
no  reasonable  approximation  to  the  response  envelopes  could  be  made.  This 
was  corrected  by  reducing  the  time  step  to  .005  sec  for  Case  IV.  Figures 
3,  6,  9,  and  12  present  the  response  envelope  plots  for  the  four  test 
cases. 

Integrated  control  effort  time  histories  for  each  simulation  were 
also  plotted.  As  expected,  these  curves  asymptotically  approach  maximum 
values  as  the  decreasing  magnitude  of  the  response  requires  reduced 
expenditure  of  control  effort.  These  curves  are  smooth  with  no  real 
scatter  in  the  data,  so  no  data  are  shown.  The  control  response  curves 
for  the  test  cases  are  presented  in  Figures  4,  7,  10,  and  13. 

By  combining  data  from  these  two  curves,  curves  showing  the  control 
configuration  performances  were  formed.  This  was  done  by  taking  sections 
of  the  response  envelope  and  control  response  plots  at  a  given  time  as 
illustrated  on  Figures  3  and  4.  These  plots  contain  families  of  curves 
with  ttQ  as  the  parameter.  Using  the  sections,  the  deflection  time  at 
time  t  may  be  plotted  versus  the  integrated  control  force  at  time  t  for 
each  value  of  cr^  used.  In  the  absence  of  prespecified  performance 
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70  degrees  and  Controller  3  had  no  relative  angle  less  than  60  degrees. 
Though  the  spillover  terms  were  significantly  reduced  and  the  systems 
were  stable. 

The  remaining  test.  Case  II,  used  no  transformation.  While  the 
individual  controllers  were  designed  for  asymptotic  stability,  no  effort 
was  made  to  insure  total  system  stability.  Case  II  did  form  a  stable 
system,  but  this  cannot  be  attributed  solely  to  its  favorable  modal 
groupings.  Case  III,  using  very  incompatible  groupings,  exhibited  total 
system  stability  before  the  transformation  was  performed. 

Though  each  test  case  had  the  possibility  of  major  problems,  all 
configurations  generated  stable  controllers  that  successfully  reduced 
system  response.  Though  triangularization  of  the  system  was  not 
achieved  in  two  of  the  cases  where  it  was  attempted,  the  magnitude  of 
the  affected  terms  were  greatly  reduced  and  successful  control  systems 
were  designed.  Simulated  response  performance  curves  were  found  for 
the  four  test  cases.  The  performance  curves  can  be  compared  two  ways. 

A  horizontal  section  taken  at  the  required  level  of  response  allows 
comparison  of  the  system  efficiencies  by  comparison  of  the  expended 
control  effort.  Taking  a  vertical  section  at  the  allowable  control 
expenditure  allows  efficiency  comparison  by  comparison  of  the  achieved 
responses.  The  first  of  these  methods  will  be  used  in  the  discussion. 
This  investigation  compares  two  groups  of  these  curves. 

The  first  group  consists  of  Cases  I  through  III.  This  forms  a  com¬ 
parison  of  three  control  configurations  applied  to  the  nominal  design 
to  explore  the  relationship  between  the  compatibility  of  the  modal 
groupings  and  triangularization.  Figure  15  presents  the  performance 
curves  for  these  three  cases.  Three  results  can  be  seen  in  Figure  15. 
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The  first  is  that  there  is  little  differentiation  between  the  configura¬ 
tions  over  a  response  range  covering  roughly  the  upper  half  of  the 
uncontrolled  response.  In  fact,  two  configurations  are  virtually  identi¬ 
cal  over  this  range.  If  the  difference  between  Case  I  and  Cases  II  and 
III  is  shown  to  be  significant  on  more  complex  models,  then  the  optimum 
control  configuration  is  shown  to  depend  on  the  required  level  of  response 
reduction.  .  .  - 

The  second  result  is  that  triangularization  is  not  desirable  if  the 
modal  grouping  is  very  compatible.  This  may  be  due  to  the  manner  in  which 
the  transformation  attempts  to  orthogonal ize  the  controllers.  It  is 
possible  that  the  transformation  may  have  predetermined  orientations  for 
the  controllers.  In  attempting  to  force  this  alignment,  the  natural 
orthogonality  of  the  favorable  grouping  is  ignored  and  more  control 
effort  is  required.  This  conjecture,  though  presently  unsupported,  could 
explain  the  worsening  of  a  favorable  grouping  by  the  transformation. 

Finally,  Figure  15  shows  the  favorable  grouping  without  triangulari¬ 
zation  to  be  slightly  more  efficient  than  the  unfavorable  grouping.  The 
inability  of  the  transformation  to  completely  triangularize  Case  II  could 
account  for  this  difference.  Successful  triangularization  would  likely 
close  the  gap  to  some  degree.  It  should  be  noted  that  this  result  depends 
on  the  quality  of  compatible  groupings  that  can  be  achieved  for  a  given 
structure.  The  deflection  design  demonstrates  that  even  the  most  favorable 
groupings  can  be  quite  incompatible. 

The  second  investigation  is  intended  to  demonstrate  comparison  of  the 
"controllability"  of  two  structural  designs.  As  implied  earlier,  this 
comparison  requires  the  application  of  the  optimum  configuration  for  each 
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structure.  For  the  investigation.  Cases  II  and  IV  are  used.  The  previous 
results  show  that  Case  II  is  the  most  efficient  configuration  of  those 
tested.  Because  the  most  favorable  groupings  of  the  deflection  designed 
structure  are  so  incompatible,  the  transformation  was  used  in  Case  IV. 

While  these  configurations  may  well  not  be  optimums  for  the  structures 
compared,  they  are  the  best  available  within  the  restrictions  of  this 
study.  There  are  several  points  of  interest  in  Figure  16  which  compares 
the  performance  curves  for  Cases  II  and  IV.  Two  surprises  are  the 
linearity  of  the  deflection  design  curve  and  the  substantial  reduction 
in  the  uncontrolled  response  for  this  case.  The  linearity  of  the  curve 
is  presently  unexplained.  The  reduction  of  the  uncontrolled  response 
on  the  deflection  design  shows  the  success  of  the  design  method.  This 
comparison  shows  a  more  clear  distinction  between  low  and  high  levels 
of  response  reduction.  The  greatest  possible  benefit  of  the  frequency 
redesign  is  found  in  the  case  where  the  required  response  at  0.9  seconds 
is  above  or  only  slightly  below  its  uncontrolled  response.  In  this  case, 
the  deflection  design  would  be  far  superior  due  to  its  reduced  uncontrolled 
response. 
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FIGURE  16:  COMPARISON  OF  PERFORMANCE  CURVES,  CASES  II  AND  IV 
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VI.  CONCLUSIONS 


This  investigation  compared  the  efficiency  of  various  control  system/ 
structural  design  configurations.  Emphasis  was  placed  on  exploration  of 
the  conditions  governing  the  loss  of  system  performance  in  transforming 
the  system  to  a  block  traingular  form.  As  expected,  the  compatibility  of 
the  modal  groupings  proved  to  be  a  significant  factor.  Comparison  of  con¬ 
troller  performance  curves  for  the  test  cases  considered  showed  that  the 
transformation  caused  significant  loss  of  system  performance  when  very 
compatible  modal  groupings  were  used.  When  an  extremely  incompatible 
grouping  was  used,  however,  the  almost  triangularized  system  performance 
compared  favorably  with  the  case  representing  an  untriangularized 
favorable  grouping.  The  triangularization  transformation  failed  due  to 
the  inadequate  number  of  sensors  and  actuators  used.  Though  the  trans¬ 
formation  failed  to  triangularize  two  test  configurations  and  was  not 
even  used  on  one  other,  all  four  configurations  generated  controllers  with 
total  system  stability. 

Examination  of  the  relative  angles  in  the  modal  groupings  of  the 
frequency  designed  structure  reveals  that  even  the  most  favorable  group¬ 
ings  for  this  design  are  very  incompatible.  The  compatibility  of  the 
modes  could  be  improved  by  relocation  of  the  actuators  and/or  sensors. 
Additional  study  is  clearly  needed  to  determine  the  degree  of  compati¬ 
bility  or  incompatibility  required  for  the  present  assessment  of  the 
three  control  configurations  to  be  valid.  Also,  a  validated  characteristic 
measure  of  grouping  compatibility  is  needed. 
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VII.  RECOMMENDATIONS 


The  results  of  this  study  have  indicated  several  areas  that  warrant 
further  exploration.  The  study  showed  on  the  CSDLI  that  spillover  elimi¬ 
nation  can  be  detrimental  or  of  benefit  to  system  performance  depending 
on  the  compatibility  of  the  modal  groupings.  The  generality  of  this 
result  should,  at  a  minimum,  be  verified  by  application  to  the  CSDLII 
model.  Further,  testing  on  more  complex  models  could  also  show  more  con¬ 
clusively  the  presence  or  lack  of  differentiation  between  the  favorable 
grouping,  no  transformation  and  the  unfavorable  grouping  with  transfor¬ 
mation  cases. 

Another  area  of  interest  is  the  degree  of  incompatibility  in  the 
modal  groupings.  A  slightly  less  incompatible  unfavorable  grouping  might 
allow  a  more  complete  elimination  of  spillover  by  the  transformation. 

This  could  eliminate  the  slight  advantage  that  the  favorable  grouping 
with  spillover  seems  to  have  at  this  point.  One  way  to  improve  the 
"worst  case"  grouping  is  to  use  additional  sensors  and  actuators  on  all 
the  test  cases.  As  has  been  noted,  the  six  pair  used  in  this  test  do 
not  guarantee  complete  elimination  of  spillover  by  the  transformation. 
This  lack  was  evidenced  in  the  unfavorable  grouping  by  the  presence  of 
zero  and  indefinite  singular  values.  The  extra  observability  and  con¬ 
trollability  would  likely  be  of  the  most  benefit  to  the  unfavorable 
grouping. 

In  addition  to  these  ideas,  the  computer  program  has  a  few  areas  of 
need  that  should  be  mentioned  here.  At  present,  the  variable  lists  in 
the  COMMON  statements  of  the  program's  subroutines  are  not  identical  to 
the  lists  in  the  main  program.  This  is  believed  to  be  the  probable  cause 
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of  two  significant  problems.  Matrices  formed  by  modifications  in  the 
program  are  not  printed  properly  by  the  subroutines  PRNT,  PRNTXL,  and 
PRNTSM.  This  has  resulted  in  the  use  of  extra  PRINT  statements  in  the 
program  that  are  not  conducive  to  changes  in  the  order  of  the  control 
model.  Also,  subroutine  MMUL  was  found  to  print  erroneous  results  with 
no  warning.  This  was  corrected  by  use  of  the  IMSL  routine  VMULFF  in 
the  modifications.  It  is  believed  that  correction  of  the  COMMON  state¬ 
ments  will  correct  these  problems  and  allow  use  of  these  subroutines. 

These  routines  are  very  important  because  their  use  allows  the  system 
order  to  be  changed  by  altering  only  the  matrix  dimension  statements. 

A  comment  section  immediately  following  the  program  initializations  could 
also  be  helpful  in  adapting  the  program  to  a  different  model.  The  needed 
comments  would  show  the  matrix  dimensions  required  in  terms  of  system 
order,  the  number  of  modes  in  each  controller,  and  the  number  of  sensors 
and  actuators.  Proper  dimensions  not  only  saves  core  memory,  but  allows 
standardized  use  of  the  VMULF_  routine  which  require  input  of  the  actual 
row  dimension  of  the  matrices.  Also,  the  needed  value  of  each  initialized 
variable  should  be  described.  A  program  of  this  size  does  not  warrant  a 
User's  Manual,  but  this  sort  of  description  section  could  be  most  helpful  - 
especially  as  the  program  is  successively  modified. 
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0.0 

.3536 

-.3536 

0.0 

0.0 

8 

0.0 

0.0 

.6124 

-.6124 

0.0 

0.0 

9 

0.0 

0.0 

.7071 

.7071 

0.0 

0.0 

10 

0.0 

0.0 

0.0 

0.0 

-.7071 

.7071 

11 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

12 

0.0 

0.0 

0.0 

0.0 

.7071 

.7071 
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PROGRAM  ACOSSCF ( INPUT , OUTPUT , T APEB » T APE6 ) 

‘  PFN:  ACOSSCF 

THIS  PROGRAM  GENERATES  A  LOWER  TRIANGULAR  TRANSFORMATION 

THE  THREE  CONTROLLER  SOLUTION  MAY  INCLUDE  RESIDUALS, 

REAL  A1<12,12),A2<12,12),A3(12,12>,A4<12,12) 

REAL  B1<12,12),B2<12,12),B3(12,12) ,B4<12,12) 

REAL  Cl <12, 12) ,C2<12,12) ,C3(12,12) ,C4<12,12) 

REAL  CTCC1<12, 12) ,CTCC2< 12,12 ) ,CTCC3( 12,12 ) ,CTCC4< 12, 12) 
REAL  SAT(12,12) ,SAT2< 12,12) ,SAT3(12,12) ,SAT4(12,12) 

REAL  AKC( 12,12) , ACT < 12, 12 ) , BCG ( 12 , 12 ) ,KCC(12,12) 

REAL  P<12,12) ,S<12,12) 

REAL  QA1(12,12) ,QA2( 12, 12) ,QA3( 12, 12) ,QA4(12,12) 

REAL  XHAT<2,48) ,PETA<21 ,48) ,CV(2,21 ) ,XL0S(2) 

REAL  ACG1<12,12),ACG2<12,12),ACG3(12,12) ,ACG4(12,12) 

REAL  ABGl <12,12), ABG2  <12,12), ABG3  <12,12), ABG4  <12,12) 

REAL  GAIN1<12,12) ,GAIN2< 12,12) , GAIN3 < 12 , 12 ) , GAIN4 ( 12 , 12 ) 
REAL  KT1<12,12) ,KT2<12,12) ,KT3<12,12) ,KT4<12,12) 

REAL  KOB1<12,12),KOB2<12,12),KOB3<12,12) ,K0B4<12,12) 

REAL  GAMMA1<12,12) ,GAMMA2<12,12),GAMMA3<12,12) 

REAL  T2< 12,12) ,T3< 12,12) ,T4< 12,12 ) ,BB< 17) 

REAL  TRT< 12,12) ,TEN< 12,12) ,CT< 12, 12) ,U< 12, 12) 

REAL  RK<12,12) ,RK1<12,12),RK2<12,12) ,RK3<12,12) 

REAL  RG2<12,12) , RG3 < 12, 12) ,RG4 < 12 , 12 > 

REAL  MAJM<48,48),D<17),W<17),T0L,DT,X0<48) 

REAL  ZETA,CL0S<2,12) ,SING<12) ,XTR < 12 , 12 ) ,X1 < 48) 

REAL  EAT<48,48) ,EAT2< 48,48), WORK <48,48 ),STOR< 12, 12) 

REAL  PHIA<12,12),PHIS<12,12),M0DE<2,12),INIT<4,12) 

REAL  PXHT<24,48) ,GPART<6,24) ,UHAT<6,48) , SCALED 
REAL  U1 <  6 ) , UBAR  <  A ) , TEFF  <  6 ) , DUMl < 1 2 ) , DUM2  <  A ) , DUM3  <  A ) 

REAL  PHI<12,12),DMAT<12,12),CMAT<12,12) ,NEWCON 
REAL  KFIRST<12,12) ,KREAL<12,12) 

REAL  PHI1< 12,6) ,PHI2< 12,6) ,PHI3< 12,6) ,PHI4< 12,6) 

REAL  PHINEW<12,12) ,GNIT<4,12) ,DUM4< 12, 12) ,DUM6<6) 

REAL  CDOT,MAGI ,CANG,PHIT< 12, 12) ,DUM5<24) 

INTEGER  N , N2 , NCI , NC2 , NC3 , NCI 2 , NC22 , NC32 , NR , NR2 , NMODE 
INTEGER  IC1<12) ,IC2<12),IC3<12),I, J,K,L,M,KK,LL,MM 
I NTEGER  DEC , Q , NACT , NSEN , IR < 12 ) , I ER , SKIP , NCOL , NCOL 1 , NM0D2 
INTEGER  NDA,NDIM,NDA1,NDIM1,ZZ,E2,E3,E4,P1,P2,P3,AA 
INTEGER  INCOND 
COMPLEX  Z<24),W1<12) 

COMMON/MAINA/NDA,NDAl 
COMMON/MAINB/NCOL , NCOL 1 
COMMON/MAINl/NDIM,NDIMl, TEN, X< 2304) 

C0MM0N/MAIN2/ST0R 
C0MM0N/MAIN3/XTR 
C0MN0N/SAUE/T<100) ,TS<100) 

COMMON/INOUT/KOUT,TAPE 

COMMON/NUM/ICl , IC2 , IC3 , IR ,NC1 , NC2 , NC3 , NR 


INITIALIZATIONS 
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PROGRAM  ACOSSCF  < INPUT  r  OUTPUT , TAPE8  r  T APEA ) 

*  pfn:  acosscf 

THIS  PROGRAM  GENERATES  A  LOWER  TRIANGULAR  TRANSFORMATION 

THE  THREE  CONTROLLER  SOLUTION  MAY  INCLUDE  RESIDUALS. 

REAL  Al<12,12),A2(12,12),A3<12rl2),A4<12,12) 

REAL  Bl<12,12) »B2( 12,12 ),D3( 12,12) ,64(12,12) 

REAL  Cl<12,12) ,C2(12,12) ,C3(12,12) ,C4(12,12) 

REAL  CTCC1(12,12) ,CTCC2( 12,12 ),CTCC3( 12,12 ) ,CTCC4< 12,12) 
REAL.  SAT<12,12) ,SAT2< 12,12), SAT3< 12,12), SAT4< 12, 12) 

REAL  AKC(12,12) ,ACT(12,12) ,BCG<12,12) ,KCC(12,12) 

REAL  P<12,12),S<12,12) 

REAL  QA1(12,12) ,QA2(12,12) ,QA3( 12, 12) ,QA4( 12, 12) 

REAL  XHAT<2,48) ,PETA<21 ,48) ,CV(2,21 ) ,XL0S(2) 

REAL  ACG1<12,12),ACG2<12,12),ACG3(12,12),ACG4(12,12) 

REAL  ABG1<12,12),ABG2(12,12),ABG3(12,12),ABG4(12,12) 

REAL  GAIN1(12,12),GAIN2(12,12),GAIN3(12,12),GAIN4(12,12) 
REAL  KT1<12,12) ,KT2< 12,12) ,KT3( 12,12 ) ,KT4( 12, 12) 

REAL  K0B1<12,12) ,K0B2<12,12) ,KOB3<12,12) ,KOB4(12,12) 

REAL  GAMMAl (12,12), GAMMA2 (12,12), GAMMA3 (12,12) 

REAL  T2(12,12),T3(12,12),T4(12,12),BB(17) 

REAL  TRT(12,12),TEN(12,12),CT(12,12) ,V(12,12) 

REAL  RK(12,12) ,RK1 ( 12, 12) ,RK2( 12, 12) ,RK3(12,12) 

REAL  RG2(12,12) ,RG3(12,12) ,RG4(12,12) 

REAL  MAJM(48,48),D<17),U(17),T0L,DT,X0(48) 

REAL  2ETA,CL0S(2,12),SING(12),XTR(12,12),X1(48) 

REAL  EAT (48, 48) ,EAT2(48,4a) ,W0RK(48,48) ,STOR( 12, 12) 

REAL  PHIA(12,12) ,PHIS( 12, 12) ,MODE( 2, 12) , INIT (4 , 12) 

REAL  PXHT(24,4a) ,GPART(6,24) ,UHAT(6,48) , SCALED 
REAL  Ul(6) ,UDAR(6) , TEFF ( 6) , DUMl ( 12 ) , DUM2 ( 6 ) , DUM3 ( 6 ) 

REAL  PHI(12,12),DMAT(12,12),CMAT(12,12) ,NEWCON 
REAL  KFIRST(12,12) ,KREAL(12,12) 

REAL  PHI 1( 12,6) ,PHI2( 12, 6),PHI3( 12,6) ,PHI4( 12,6) 

REAL  PHI NEW (12,12),GNIT(4,12), DUM4 (12,12), DUM6 ( 6 ) 

REAL  CDOT , MAGI , CANG , PHI T ( 1 2 , 12 ) , DUM5 ( 24 ) 

I NTEGER  N , N2 , NCI , NC2 , NC3 , NCI 2 , NC22 , NC32 , NR , NR2 , NMODE 
INTEGER  IC1(12) ,IC2(12) ,IC3(12),I, J,K,L,M,KK,LL,MM 
INTEGER  DEC , 0 , NACT , NSEN , IR ( 12) , lER , SKIP , NCOL , NCOLl ,NM0D2 
INTEGER  NDA,NDIM,NDA1,NDIM1,ZZ,E2,E3,E4,P1,P2,P3,AA 
INTEGER  INCOND 
COMPLEX  Z(24),W1(12) 

COMMON/MAINA/NDA,NDAl 
COMMON/MAINB/NCOL , NCOLl 
C0MM0N/MAIN1/NDIM,NDIM1,TEN,X(2304) 

C0MM0N/MAIN2/ST0R 
C0MM0N/MAIN3/XTR 
C0MM0N/SAUE/T(100) ,TS(100) 

COMMON/INOUT/KOUT,TAPE 

C0MM0N/NUM/IC1,IC2,IC3,IR,NC1,NC2,NC3,NR 
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INITIALIZATIONS 


^U  ^L  «b 

9  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  f|(  ^  ^ 

READ (8,#)  SCALEB 

^u 

^U  ^. 

^A  M  •  ^n  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  •^ 


NDIM  =  12 
NDIMl  =  13 
NCOL  =  12 
NCOLl  =  13 
NM0ri2  =  24 
NDA  ==  48 
NDAl  =  49 
KOUT  =  6 
TAPE  =  9 
0  s  0 
lER  =  0 
2Z  =  0 
AA  «  0 


PRINT' <////) ' 

tj^  p%  9  ^1  y  ^1  ^u  ^u  kW  ^p 

JL  1^  1  ^A  y  ^A  4A  ^n  ^A ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A ^A ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  4A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  *^A 


PRINT*,' 

***** 

***** 

PRINT*, ' 

***** 

34L0WER-RESIIIUAL 

***** 

PRINT*,' 

***** 

BLOCK 

***** 

PRINT*, ' 

***** 

C  s  n  L  11 

***** 

PRINT*, ' 

***** 

***** 

PRINT*, ' 

***** 

CONTROL  FORCE  MODIFIED 

^  ^  ^  y 

^  ^  ^  ^ 

PRINT*, ' 

***** 

]}( 3|(  )|C  )|C  3|(  ^ 

p%  V  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^u  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  y«  ^  ^  ^  ^  ^  ^  ^ 

0  I  y  ^A  ^A  ^A  ^n  ^A  ^n  ^A  ^A  ^n  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A  ^A 

PRINT' <//)  ' 

PRINT*,'  THIS  PROGRAM  GENERATES  A  SOLUTION', 

*  '  USING  A  LOWER  TRIANGULAR  TRANSFORMATION  ' 

PRINT'<////)' 


INITIAL  SELECTION  FOR  THREE  OR  FOUR  CONTROLLERS 


oooo  ooo  o  o  o  noo 


PHI  MiMRICES  rtNIi  CONTROLLER  ENTRIES 


PRINT' (///) ' 

IF  <DEC.EQ.3)  THEN 

PRINT*,'  ENTER  NCI , NC2, NC3,NR, NACT , NSEN , ZETA  >' 
ELSE 

PRINT*,'  ENTER  NCI , NC2 , NC3 ,NC4 ,NACT , NSEN , ZETA  >' 
END  IF 

READ ( 8 , * )  NCI , NC2 , NC3 , NR , NACT , NSEN , ZET A 
PRINT* , NC 1 , NC2 , NC3 , NR , NACT , NSEN , ZETA 
PRINT*,'  ' 

N  =  NCI  +  NC2  +  NC3  +  NR 
READ  PHI  MATRIX  BY  COLUMNS 
PRINT' (/) ' 

DO  1  J=1,N 

READ(8,»)(PHI(I,J>,I=1,N) 

PRINT*, 'MODE  SHAPE  ',J,'  IS>' 

PRINT*,'  ' ,<PHI<I,J),I=1,N) 

1  CONTINUE 
PRINT' <//) ' 

READ  THE  D  MATRIX  BY  ROWS 

PRINT*, 'THE  D  MATRIX  IS' 

PRINT' </) ' 

DO  2  1  =  1, N 

READ<8,*)  <DMAT<I, J) , J=1 ,NACT) 

PRINT' <1X,12F10.4) ' , <DMAT( I , J) , J=l ,NACT) 

2  CONTINUE 
PRINT' <//)' 

READ  CMAT  TRANSPOSED  BY  ROWS 

PRINT*, 'THE  C  MATRIX  TRANSPOSED  IS' 

PRINT' </) ' 

DO  3  1  =  1, N 

READ<8,*)  <CMAT<I,J) ,J=1,NSEN> 

PRINT' <1X,12F10.4) ' , <CMAT( I , J) , J=l ,NSEN) 

3  CONTINUE 
PRINT' <//) ' 


NOW  COMPUTE  PHI A  AND  PHIS  <PHICTR3*D  AND  C*PHI) 

CALL  VMULFM < PHI , DMAT ,N ,N ,NACT , NCOL , NCOL , PHI A , NCOL , lER) 
PRINT*, 'THE  PHIA  MATRIX  IS' 

CALL  PRNT(PHIA,N,NACT) 

CALL  T ANGL ( PHI A , N , NACT , DUMl , DUM4 ) 

PRINT' (//)  ' 

CALL  VMULFM < PHI , CMAT , N , N , NSEN , NCOL , NCOL , PHIS , NCOL , lER > 
PRINT*, 'THE  PHIS  MATRIX  IS' 

PRINT' </) ' 

CALL  PRNT<PHIS,N,NSEN) 

PRINT' <//) ' 


OMEGAS 
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C  INITIAL  CONDITIONS 

C  THESE  ARE  READ  IN  ONLY  ONCE  FOR  EACH  JOB, 

C  REGARDLESS  OF  THE  NUMBER  OF  CONSECUTIVE  RUNS 
C 

C  REMEMBER  —  EVEN  THOUGH  RESIDUAL  MODES  DO  NOT  HAVE  E  AND 
C  E  DOT  TERMS,  THIS  READ  STATEMENT  WILL  BE  LOOKING  FOR  THEM, 

C  SO  INPUT  E  AND  E  DOT  TERMS  FOR  THE  RESIDUALS  ALSO, 

C  DON'T  WORRY  —  SUBROUTINE  FORMXO  WILL  FILTER  THEM  OUT  LATER. 

C 

C  INCOND  =  0  IF  INITIAL  CONDITIONS  ARE  IN  PHYSICAL  VARIABLES 
C  INCOND  =  1  IF  INITIAL  CONDITIONS  ARE  IN  MODAL  COORDINATES 

INCOND  =  1 

J*\ 

^p 

c 

IF(AA.EQ.O)  THEN 

PRINT*,'  ENTER  THE  INITIAL  CONDITIONS  FOR  ',N,'  MODES  ' 

IF  (INCOND. EQ.O)  THEN 

PRINT*,'  CALL  THE  PHYSICAL  VARIABLE  'G*  ' 

PRINT*, 'THEN  ROWl=G,  R0W2=G  DOT,  ROWS^^'E,  ROW4=E  DOT' 

DO  80  1=1,4 

READ(8,»)  (GNITd,  J)  ,  J=1,N) 

PRINT' < IX, 12F10.8) ' , <GNIT< 1 , J) , J=1 ,N) 

80  CONTINUE 

CALL  TFR<PHIT,PHI,N,N,1,2) 

CALL  LINV3F<PHIT,DUM1,1,N,N,-1,D2,DUM5,IER) 

CALL  VMULFF<GNIT,PHIT,4,N,N,4,N,INIT,4,1ER> 

PRINT' </) ' 

PRINT*,'  INITIAL  CONDITIONS  IN  MODAL  COORDINATES  ARE  ' 

DO  81  1=1,4 

PRINT' < 1X,12F10,7) ' , ( INIT< I , J> , J=1 ,N) 

81  CONTINUE 
ELSE 

PRINT*,'  R0W1=ETA,  R0W2=ETA  DOT,  R0W3=E,  R0W4=E  DOT  ' 
PRINT*, 'NOTE  ... .INPUT  IS  IN  MODAL  COORDINATES!' 

DO  ?0  1=1,4 

READ<8,*)  <INIT(I,J),J=1,N) 

90  CONTINUE 

CALL  VMULFP(INIT,PHI,4,N,N,4,N,GNIT,4,IER) 

PRINT' <///) ' 

PRINT*, 'INITIAL  CONDITIONS  IN  PHYSICAL  VARIABLES  ARE' 

DO  94  1=1,4 

94  PRINT' <1X,12F10,8) ' , <GNIT< I , J) , J=1 ,N) 

END  IF 

PRINT' <//) ' 

PRINT*,'  THE  INITIAL  CONDITIONS, (XO  AT  T=0),  ARE  ' 

CALL  FORMXO(XO,INIT,DEC) 

CALL  PRNT(X0,M,1) 

AA  =  1 
END  IF 

100  CONI INUE 

PRINT*,'  TO  PRINT  ALL  OF  THE  MATRICES  ENTER  1,  ELSE  ENTER  0 
READ(8,*>  Q 


ooo  o o o n n  nnnnnnnn 


READ  IN  THE  WEIGHTING  MATRIX  DIAGONAL  OALUES 


note:  the  input  values  are  for  the  purpose  of  SETTING 

ONLY  THE  RELATIVE  MAGNITUDE  OF  WEIGHTING  OF  THE  MODES. 
THESE  ARE  THEN  SCALEDD  BY  THE  VARIABLE  ’SCALEB*. 


PRINT*,'  ENTER  THE  DIAGONAL  VALUES  OF  QOB  ' 

PRINT*,'  ' 

PRINT*,'  note:  these  VALUES  WILL  BE  SCALED  BY  THE  VALUE  ' 
PRINT*,'  OF  'SCALEB'  SET  IN  THE  PROGRAM  INITIALIZATIONS  ' 
READ<8,*>  (BB<I),I=1,N) 

PRINT*,'  ' 

PRINT*, <BB(I), 1=1, N) 

DO  106  1=1, N 
106  BB<I>=BB(I)*SCALEB 
PRINT*,'  ' 

PRINT*, 'WEIGHTING  MATRIX  IS  MULT  BY  SCALE  =  ',SCALEB 
PRINT' <///) ' 


FORMING  THE  A,B,C  AND  WEIGHTING  MATRICES 


CALL  F0RMA<A1,D,W,NC1,NC12,IC1> 

CALL  F0RMB(B1 ,PHIA,NC1 ,NC12,NACT, ICl ) 
CALL  F0RMC<C1,PHIS,NC1,NC12,NSEN,ICI) 
CALL  F0RMPHI<PHI1,PHI,N,NC1,IC1) 

CALL  FORMA ( A2, D, W, NC2, NC22, IC2> 

CALL  FORMB ( B2 , PHI A , NC2 , NC22 , NACT , IC2 ) 
CALL  F0RMC(C2,PHIS,NC2,NC22,NSEN, IC2> 
CALL  F0RMPHI(PHI2,PHI,N,NC2,IC2> 

CALL  F0RMA(A3,D,W,NC3,NC32,IC3) 

CALL  FORMB <  B3 , PHI A , NC3 , NC32 , NACT , IC3 ) 
CALL  F0RMC(C3,PHIS,NC3,NC32,NSEN,IC3> 
CALL  F0RMPHI(PHI3,PHI,N,NC3,IC3) 

CALL  FORMA <A4,D,W, NR, NR2,IR) 

CALL  F ORMB  <  B4 , PHI A , NR , NR2 , NACT , I R ) 
CALL  FORMC  < C4 , PHI S , NR , NR2 , NSEN , I R ) 
CALL  F0RMPHI<PHI'1,PHI,N,NR,IR) 

CALL  F0RMQ1(0A1,BB,NC1,IC1) 

CALL  F0RMQ1(0A2,BB,NC2,IC2) 

CALL  F0RMQl(aA3,BB,NC3,IC3) 

IF  <DEC.Ea.4)  THEN 
CALL  F0RMQ1(GA4,6B,NR,IR> 

END  IF 


PRINTING  THE  A,B,C  AND  WEIGHTING  MATRICES 
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EAT2  IS  NOW  THE  SOLUTION  TO  E  TO  THE  AT 

WE  NOW  PROPAGATE  THE  E  TO  THE  AT  MATRIX  IN  HALF  MINUTE  STEPS 
FIRST  GENERATE  THE  CCi;C2:C33  MATRIX 


DO  260  I--=1,NSEN 
no  260  J=1,NC1 
260  V(I, J)=C1(I,J) 
no  270  I=^U,NSEN 
no  270  J=1,NC2 
270  V<I, J+NC1)=C2<I, J) 
no  280  I=1,NSEN 
no  280  J=1,NC3 
280  V<I, J+NC1+NC2)=C3(I,J) 
NM0nE=NCl+NC2+NC3 
IF  (NR.GT.O)  THEN 
no  290  I=1,NSEN 
no  290  J=1,NR 
290  U<I,J4NM0IiE)=C4(I,  J) 

NMOnE=NMOnE+NR 

ENHIF 

PRINT*,'  THE  C  PARTITION  MATRIX  IS  ' 
CALL  PRNT<V,NSEN,NMOnE> 

NOW  GENERATE  THE  CPHIl :PHI2:PHI3D  MATRIX 

no  282  1=1, N 
no  282  J=1,NC1 
282  PHINEW(I,J)=PHI1<I,J) 
no  284  1  =  1, N 
no  284  J=1,NC2 

284  PHINEWd,  J+NC1)=PHI2<I,J) 

no  286  1  =  1, N 

no  286  J=1,NC3 

286  PH1NEW<I,J+NC1+NC2)=PHI3<I, J> 
NM0nE=NCl+NC2+NC3 
IF (NR.GT.O)  THEN 
no  288  1=1, N 

no  288  J=1,NR 

288  PHINEWd,  J+NM0nE)=PHI4<  I,  J) 
NMOnE=NMOnE+NR 
ENniF 

PRINT*,'  THE  REARRANGEIi  PHI  MATRIX  IS' 
CALL  PHNT<PHINEW,N,NMOnE) 


NEXT  GENERATE  THE  CGI :G2:G3:G43  MATRIX 


non  no  oo  o  non 


NEXT  GENERATE  THE  CGI :G2:G3:G43  MATRIX 


DO  291  I=1,NSEN 
DO  291  J=1,NC12 

291  GPART<I,J)=GAIN1(I,J) 

DO  292  I=1,NSEN 

DO  292  J=1,NC22 

292  GPART<I»J+NC12)=GAIN2<I,J) 

DO  294  I=1,NSEN 

DO  294  J=1,NC32 

294  GPARTd  ,  J+NC12+NC22)=GAIN3(I,  J) 
NM0D2=NC12+NC22+NC32 
IF  (NR.GT.O)  THEN 
DO  296  I=1,NSEN 
DO  296  .J=1,NR2 

296  GPARTd  ,J+NM0D2)=GAI  N4<  I,  J) 
NM0D2=NM0D2+NR2 
ENDIF 

PRINT#,  dHE  6  PARTITION  MATRIX  IS' 

DO  66  I:=1,NSEN 

66  PRINT'  <  IX,  12F10.5)  '  ,  (GPARTd,  J)  ,  J=l,12) 
PRINT' <//) ' 

DO  67  I=1,NSEN 

67  PRINT' <1X,12F10.5) ' , < 6PART d , J> , J=13 ,24 > 
PRINT' <///) ' 


CALL  THAT  <  CLOS , NSEN , NCI , NC2 , NC3 , NR , MM , XHAT , PET A , PXHT , CV , V 
i,NMOD2,PHINEW) 

PRINT#,'  THE  SCALING  FACTOR  IS  ', SCALED 

CALL  VMULFF  <  GPART , PXHT , NSEN , NM0D2 , MM , NSEN , NM0D2 , UHAT , NSEN , I ER  J 
CALL  T I  ME  <  EAT2 ,  MM ,  DT ,  NACT ,  XO  ,  XI ,  XHAT ,  UHAT ,  EAT ,  WORK ,  liUM6 , 
JDUM2,DUM3) 

PRINT' (//) ' 


300  CONTINUE 


EIGENVALUE  ANALYSIS  SECTION 

?? 
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SUBROUT  I NE  TANGL  ( F  I E  A » N  »  NACT ,  F-MAG ,  ANGL ) 

REAL  FIEA(N,NACT> ,ANGL<N,N) ,PHAG(N) 

■real  CD0T»MAGI,CANG 
INTEGER  I,J,KfN,NACT 

CONFUTE  MAGNITUDE  OF  EACH  MECTOR 
DO  2  1=1 fN 
MAGI  =  0*0 
DO  1  J=1,NACT 

MAGI  =  MAGI  +  (FIEAd ,  J)**2) 

1  CONTINUE 
PMAG<I)  =  <MAGI)K*0,5) 

2  CONTINUE 

C  COMPUTE  DOT  PRODUCTS 
DO  5  1=1, N 
DO  A  J=I,N 
CDOT  =  0.0 
DO  3  K=1,NACT 

CDOT  =  CDOT  +  FIEA<I,K)*FIEA<J,K> 

3  CONTINUE 
ANGL<I r J)=CDOT 

4  CONTINUE 

5  CONTINUE 

C  COMPUTE  THE  ANGLES 
DO  7  J=lrN 
DO  6  I=J,N 

GANG  =  ANGL<J,I)/(PMAG<I)*PMAG( J)> 

IF<CANG.GT.1.0)THEN 
GANG  =  1.0 
END  IF 

ANGLd  ,  J)=ACOS<CANG) 

ANGL<I,J)=ANGL(I,J)  »  57.295779514 
IF<ANGL<I, J) .GT.90. )THEN 
ANGLdfJ)  =  ANGLd,J)  -  180.0 
END  IF 

6  CONTINUE 

7  CONTINUE 

C  PRINT  ANGLE  TABLE 
DO  9  J=1,N 
DO  8  I«J+1,N 
ANGL<  J,I)=ANGLd,  J) 

8  CONTINUE 

9  CONTINUE 
PRINTM/) ' 

PRINT*,'  THE  RELATIVE  ANGLES  BETWEEN  THE  MODES  ARF 

PRINT*,' 

PRINT*,'  MODE  1234 

t  7  8  9  10  11  12' 

DO  10  1=1, N 

PRINT' < IX, 14, 12(2X,1F7.2) )',!,< ANGL< I, J>,J=1,N> 

10  CONTINUE 
PRINT' </) ' 

RETURN 


nnnonnno 


SUBROUTINE  YHrtT  < CLOS ,  NSEN ,  NCI  r  NC2 r  NC3 ,  NR , MM r  XHAT  r  PETA  r  PXHT ,  CV  v 
»,NM0Il2rPHINEW) 

■  COMMON/MAINA/NDA,NIiAl 
COMMON/MAINB/NCOL 
INTEGER  L,LL,T,TT,K3 

INTEGER  I  r  J ,  K ,  KK ,  MM ,  NSEN ,  NACT ,  NMODE ,  NM0rt2 

REAL  CV<2,NC0L)  ,y<NCOL»NCOL)  ,PElA(NCOL,NliA)  ,PXHT<NM0i;i2,NriA) 
REAL  XHAT(2,NtiA>  ,CL0S<2»NC0L)  ,PHINEW(NCOL,NCOL) 
NM0DE=NC1+NC2+NC3+NR 

CALL  VMULPP < CLOS , PHINEU , 2 , NMODE , NMODE , 2 , NCOL , CU , 2 , lER ) 

PRINT*,'  THE  CLOS  *  PHINEW  PARTITION  MATRIX  IS  ' 

PRINT*,' 

DO  76  1=1,2 

76  PRINT' < IX, 12F10.5) ' , (CV( I , J> , J=1 , 12> 

PRINT' (//) ' 

GENERATE  THE  MATRIX  P<ETA)  TO  EXTRACT  ETA 
FROM  THE  FULL  STATE  VECTOR  Z 


THE  MATRIX  PXHT  EXTRACTS  XHAT  < ETAHAT : ETADOTHAT > 
FROM  THE  FULL  STATE  VECTOR  Z 

DO  60  1=1, NMODE 
DO  59  J=1,MM 
PXHT<I+NM0DE,J)=0.0 
PXHTd,  J)=0.0 

59  PETAd,  J)=0,0 

60  CONTINUE 
L=2*NC1 
LL-L+NCl 

DO  61  1=1, NCI 
PETA<I,I)=1.0 
PXHTd, I)  =  1,0 
PXHTd,I+L)  =  l  .0 
PXHTd+NCl,I+NCl)  =  l,0 
PXHTd+NCl,I+LL)  =  1.0 

61  CONTINUE 
K=NC1 
T=2*NC1 
KK=4*NC1 
K3=KK+NC2 
L=KK+2*NC2 
LL»L+NC2 
TT=T+NC2 

DO  62  I=1,NC2 
PETAd+K,I+KK)  =  1.0 
PXHTd+T,I+KK)  =  1.0 
PXHrd+T,I+L)-1.0 
PXHTd+TT,I+K3>  =  l,0 
PXHT<I+TT,I+LL)=1.0 

62  CONTINUE 
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K=K+NC2 
T=T+2*NC2 
•  KK=KK+4#NC2 
K3=KK+NC3 
L=KK+2*NC3 
LL*L+NC3 
TT=T+NC3 
DO  63  I»lfNC3 
PETrt(I+K,HKK)  =  1.0 
PXHT<I+T,I+KK)=1.0 
PXHT<I+T,I+L)=1 .0 
PXHT<I+TT,I+K3)=1.0 
PXHT<HTT,I+LL)  =  1 .0 
i  CONTINUE 

IF  (NR.EQ.O)  GOTO  65 

K=K+NC3 

T=T+2»NC3 

KK=KK+‘4*NC3 

K3=KK+NR 

TT=T+NR 

L=KK+2»NR 

ll=l;nr 

DO  64  1=1 rNR 

PETA<I+K,I+KK)=1.0 

PXHT<I+Trl+KK)=1,0 

PXHT<I+T,I+L)=1.0 

PXHT(I+TT,I+K3)=1*0 

PXHT<I+TT,I+LL)=1.0 

CONTINUE 

CONTINUE 

PRINT*,'  THE  MATRIX  P<ETA)  IS  ' 

PRINT' </) ' 

DO  66  I=1,NM0DE 

66  PRINT' <1X,12F10.5)' ,<PETA(I,J),J=1,12> 
PRINT' <//) ' 

DO  67  I=1,NM0DE 

67  PRINT'<1X,12F10.5)' ,<PETA<I,J),J=13,24) 
PRINT' <//) ' 

DO  68  I=1,NM0DE 

68  PRINT' <1X,12F10.5)' ,<PETA<I,J>,J=25,36) 
PRINT' (//)' 

DO  69  I=1,NM0DE 

69  PRINT' <1X,12F10.5) ' , <PETA< I , J) , J=37,48) 
PRINT' <////) ' 

PRINT*, 'THE  MATRIX  P<XHAT)  IS' 


n  o  o  n  n 


PRINT*, 'THE  MATRIX  P<XHAT)  IS' 

.  PRINT'</>' 

DO  71  I=l,NM01i2 

71  PRINT' (1X,12F10.5) ' , <PXHT( I , j> , J=1 , 12) 

PRINT' <//) ' 

DO  72  I«1,NM0D2 

72  PRINT' (IX, 12F 10.5) ' , < PXHT < I , J) , J=13 ,24 ) 

PRINT' <//) ' 

DO  73  I-l,NM0Ii2 

73  PRINT' <1X,12F10.5) ' , (PXHT < I , J) , J=2S , 36 > 

PRINT' (//) ' 

DO  74  I»1,NM0D2 

74  PRINT' (1X,12F10.5) ' , ( PXHT ( I , J) , J=37 , 48 ) 

PRINT' (////) ' 

CALL,  VMULFF  ( C  V ,  PET  A ,  2 ,  NMODE  ,  MM  ,  2 ,  NCOL ,  XHAT ,  2 , 1 ER ) 
PRINT*, 'XHAT' 

DO  75  1=1,48 

PRINT' (1X,2F10.6) ' , ( XHAT ( J, I ) , J=1 , 2 ) 

75  CONTINUE 

XHAT  =  CLOS  *  PHINEW  *  P(ETA) 

UHAT  =  CG1:G2:63}G4]  *  P(XHAT) 


RETURN 


SUBROUTINE  T I ME ( EAT2 , MM » DT , NMCT , XO , X 1 » XHAT , UHMT , E AT , WORK 
8UEiAR,TEFF,Ul) 

‘  COMMON/MAINA/NDA,NriAl 
COMMON/INOUT/KOUT 

REAL  A,ZrE;AT2(NDA,NliA)  ,EAT(NDA,NriA)  , WORK ( NDA , NDA > 

REAL  DT,XO(NDA) ,X1<NDA) ,XHAT(2,NDA) ,XL0S(2) ,NEWCON 
REAL  UBAR<NACT) ,TEFF<NACT) »UHAT(NACT,NDA) ,U1(NACT) 
INTEGER  I, J,K»M,KK,MM,NACT 
KK  «=  0 

BO  5  I=1,NACT 
5  TEFF<I)=0.0 
NEUCON=0.0 
A  *  0. 

Z  =  DT  »  1. 

PRINT#, ' 

I  INTEGRATED  EFFORT  OF  ACTUATOR  AT  TIME  T' 

PRINT*, ' 

*1  2  3  4  5  T 

PRINT#.'  TIME  XLOS  YLOS  RADIUS 


I - 

DO  10  1=1, MM 
DO  10  J=1,MM 
•  WORK<I,J)  =  EAT2<I,J) 

CALL  UMULFF(XHAT,X0,2,48,1 ,2,48,XL0S,2,IER> 

CALL  VMULFF(UHAT,X0,NACT,NDA,1,NACT,NDA,UBAR,NACT,IER) 

DO  6  I=1,NACT 
6  U1<I)»UBAR<I) 

RAD  =<<XL0S<1)##2)  +  <XL0S<2)##2> >##0.5 
WRlTE<#,15)A,XLOS<l),XLOS<2),RAri,TEFF(l),TEFF(2),TEFF(3) 
*TEFF<4) ,TEFF<5) ,NEWCON 

5  F0RMAT<2X,1F5.3,3<2X,1F10.8) ,10X,5(1X,F10.B) ,1X,1F11 .8) 

CONTINUE 
M  «  0 
CONTINUE 

IF<KK.GT.O)  THEN 

CALL  VMULFF  <  WORK , E AT2 , MM , MM , MM , NBA , NBA , EAT , NDA , lER ) 

DO  40  1=1, MM 
DO  40  J=1,MM 
EAT2<I,J)  =  EAT<I,J) 

CONTINUE 
ENDIF 
KK  =  1 

CALL  VMULFF ( E AT2 , XO , MM , MM , 1 , NDA , NDA , X 1 , NDA , I ER ) 

DO  50  1=1, MM 
XO<I)  «  X1<I> 

NEWCON  =0.0 

CALL  VMULFF  <  UHAT , XO , NACT , NDA , 1 , NACT , NDA , UBAR , NACT , I ER ) 

DO  55  1=1, NACT 

TEFF<I>«TEFF<I)  +  DT  #  ABS<UBAR(I)  +  Ul(I)>/2 
U1(I>-UBAR(I) 

NEWCON  =  NEWCON  -I-  TEFF(I) 

55  CONTINUE 


M  =  M  +  1 

IF(<M*DT) »EQ.Z)  THEN 
A  =  A  +  Z 

CALL  VMULFF ( XHAT , XO , 2 , NDA , 1 , 2 , NDA , XLOS , 2 , 1 ER  > 

RAD  =<  <XL0S(1)J|{»2)  +  <XL0S(2>*)K2)  ))t<»0.5 

WRITE<*,15)A,XL0S(1> ,XL0S<2> ,RAD,TEFF(1> ,TEFF(2) ,TEFF<3) 
STEFF<4) ,TEFF<5) ,NEUCON 
IF(A.0T.1,50)  GOTO  60 
GOTO  20 
ELSE 
GOTO  30 
END  IF 
CONTINUE 
RETURN 


SUBROUT I NE  FORMPHI ( FEE , PHI , N , NC , IC  > 

COHMON/MAINB/NCOL 

;?E AL  FEE  <  NCOL ,  NC ) , PHI  <  NCOL , NCOL  > 

INTEGER  IC(NC>,M,N,NC,I,J 

DO  2  J=1,NC 

H  =  IC<J> 

DO  2  I»1,N 
FEE(I,J>=  PHI(I,M) 

2  CONTINUE 
RETURN 
END 


SUBROUTINE  MMUL<X,Y,N1 ,N2,N3,Z) 

COMMON  /MAINB/  NCOL 

DIMENSION  X<NCOL,l) ,Y<NCOLrl) ,Z<NCOL»l> 
DO  3  J=lfN3 
DO  2  I=-1,N1 
S=0* 

DO  1  K=1,N2 

1  S*=S+X<I,K>)KY<K,J) 

2  Z<I,J)=S 

3  CONTINUE 
END 


SUBROUTINE  PRNTSM<MAT,NrM) 

COMMON/MAINB/NCOL 

REAL  MAT <2, NCOL) 

INTEGER  I,J,KrLfMfN 
PRINT*,'  ' 

DO  1  L«lfMfl2 
K  ■  L  +  11 

IF  (M-L.LT.ll)  K  »  M 
DO  2  I»1,N 

PRINT' < IX, 12F10.5) ' , <MAT<I , J) , J=L,K) 
2  CONTINUE 

PRINT'<//)' 

1  CONTINUE 

PRINT' <///)' 

RETURN 

END 


SUBROUTINE  PRNT<MAT,N,M) 
COMMON/MAINB/NCOL 
‘  REAL  MAT(NCOL,NCOL) 

INTEGER  N,I,J,K»M 
PRINT*,'  ' 

IF  (M.GT.12)  GOTO  2 
DO  1  1=1, N 

PRINT' < IX, 12F10.4) ' , < MAT < I , J) , J=1 , M > 

CONTINUE 

GOTO  10 

CONTINUE 

IF  <M.GT,24)  THEN 
CALL  PRNTXL<MAT,N,M) 

RETURN 

ENDIF 

DO  3  1=1, N 

PRINT' < IX, 12F10.4) ' , <MAT<I, J) ,J=1,12) 

CONTINUE 

PRINT' <//) ' 

DO  A  1=1, N 

PRINT' < IX, 12F10.4)' , <MAT<I , J) , J=13,M) 

CONTINUE 

PRINT'<///>' 

RETURN 

END 


SUBROUTINE  PRNTXL ( MAT ,N ,M ) 
COMMON/MAINB/NCOL 
COMMON/MAINA/NDA 
REAL  MAT(NDA,NDA) 

INTEGER  I,J,K,L,M,N 
PRINT*,'  ' 

DO  1  L=1,M,12 
K  =  L  +  11 

IF  (M-L.LT.ll)  K  =  M 
DO  2  1=1, N 

PRINT' <1X,12F10.5) ' , (MAT< I , J> , J=L,K> 

CONTINUE 

PRINT' <//>' 

CONTINUE 
PRINT' <///)' 

RETURN 

END 
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Modem  optimal  control  methods  are  used  to  develop  a  three  controller 
system.  The  control  problem  is  formulated  in  state  vector  form  using  full 
state  feedback.  ITie  control  is  implemented  on  a  modal  representation  of 
the  structural  equations  of  motion.  State  estimates  are  provided  by  a 
deterministic  observer.  A  conputer  model  is  used  to  generate  controller 
and  estimator  gains  using  steady  state  optimal  regulator  theory.  Means 
for  extracting  structural  response  and  control  effort  time  histories  are 
developed  and  related  to  the  computer  model. 

The  three  controller  system  is  used  to  control  the  CSDL  I  model 
subjected  to  physical  initial  conditions.  This  model  is  a  lumped -mass 
tetrahedral  truss  consisting  of  four  unit  masses  that  are  pin-connected 
by  massless  rods.  Two  versions  of  the  model  are  used:  the  nominal 
version  as  designed  by  the  Draper  Lab  and  a  perturbed  version  designed  to 
minimize  the  deflection  due  to  a  static  load. 

Four  structure/controller  configurations  are  tested  via  computer 
simulations.  Compatibility  of  modal  groupings  and  the  use/non-use  of 
system  triangularization  are  the  factors  varied  to  form  the  controller 
configurations.  All  configurations  use  three  controllers  each  control¬ 
ling  four  of  the  twelve  modes.  Modal  input  to  the  program  is  supplied 
by  a  NASTRAN  analysis  of  each  structure.  Time  histories  of  structural 
response  and  expended  control  effort  are  computed  for  each  test  config¬ 
uration  using  a  series  of  structural  response  weightings  to  vary  the 
control  and  estimator  gains.  Performance  curves  are  formed  for  each 
configuration  to  allow  comparison  of  system  efficiencies.  The  effect  of 
the  combination  of  modal  compatibility  and  use/non-use  of  system 
triangularization  on  system  effiency  is  examined.  The  relative 
"controllability''  of  the  two  structures  is  also  discussed. 
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